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ABSTRACT 

A correlation study has been conducted to determine the ability of cur- 
rent analytical spin prediction techniques to predict the flight motions of a 
current fighter airplane configuration during the spin entry, the developed 
spin, and the spin recovery motions. The airplane math model used aerodynam- 
ics measured on an exact replica of the flight test model using conventional 
static and forced-oscillation wind-tunnel test techniques and a recently de- 
veloped rotation-balance test apparatus capable of measuring aerodynamics 
under steady spinning conditions. An attempt was made to predict the flight 
motions measured during stall/ spin flight testing of an unpowered, radio- 
controlled model designed to be a scale, dynamically-scaled model of a 
current fighter configuration. Comparison of the predicted and measured 
flight motions showed that while the post- stall and spin entry motions were 
not we ll -predicted, the developed spinning motion (a steady flat spin) and 
the initial phases of the spin recovery motion were reasonably well predicted. 
The study clearly demonstrated the need to model the effects of spin rotation 
rate on the aerodynamics, the large influence of the pitch damping coeffeicient 
on the developed spin, and the extreme sensitivity of the predicted model mo- 
tions to control timing and phasing. 

SUMMARY 

The NASA Langley Research Center has a broad research program designed to 
advance the state-of-the-art in the area of stall/ spin technology. A major 
requirement in this area is to develop and validate reliable theoretical methods 
for the prediction and analysis of airplane stall/ spin characteristics. The 
present study was performed as the final phase of a two-part study to determine 
factors which influence the correlation between predicted and measured airplane 
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spin motions and to actually correlate analytically predicted spin motions with 
flight motions measured on a radio-controlled, dynamically- scaled drop-model of 
a current fighter airplane. This final study concentrated on the correlation 
between predicted and measured flight motions and considered the spin entry, 
the developed spin, and the spin-recovery motions. The airplane math model 
used in the study was based on the full, nonlinear equations of motion and em- 
ployed aerodynamics measured in low-speed wind-tunnel tests on an exact replica 
of the scale model used in the flight tests. These aerodynamic data were 
obtained from conventional static and forced-oscillation test techniques and 
from a recently developed rotation-balance test apparatus capable of measuring 
the effect of spin rotation rate on airplane aerodynamics. Extensive analysis 
was made of the measured aerodynamics to develop a comprehensive aerodynamic 
model of the airplane. Control time histories taken from the model flight test 
results were used in the analytical model to attempt to calculate the stall, 
spin-entry, developed spin, and recovery motions measured in the radio-con- 
trolled drop-model tests. The calculated results were correlated with the 
flight test results. 

Results of the study showed that the model developed- spin motion (a 
flat, steady spin) and the early phases of the spin-recovery motion were rea- 
sonably well predicted but that the spin-entry and complete spin-recovery mo- 
tions were not as well predicted. The high degree of correlation obtained was 
only possible when the effects of rotation rate on the model aerodynamics were 
included - that is, when rotation-balance data were used. In that regard, it 
was also noted that the phases of the motion which were not well predicted 
were the flight regions for which rotary aerodynamic characteristics had to be 
assumed. The predicted model motions were found to be very sensitive to small 
changes in the timing and phasing of the control inputs. 


The character of the 
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calculated developed spin motion was found to be sensitive to the value of the 
pitch-damping coefficient. Finally, the lack of good body attitude information 
on the flight model was found to be a severe handicap in determining the pro- 
per initial conditions to use in the analytical study since the initial at- 
titudes assumed were found to dramatically effect the spin entry motion. 


iv 


CONTENTS 

PAGE 

SUMMARY i 

INTRODUCTION 1 

TECHNICAL APPROACH - 4 

Radio Control Model Time Histories 4 

Aerodynamic Model 5 

Inertia and Mass Data 8 

Analytical Technique g 

RESULTS AND DISCUSSION 22 

CONCLUDING REMARKS 16 

REFERENCES ]_8 

LIST OF FIGURES - 19 

APPENDICES 36 

Appendix A: Equations of Motion and Associated Formulas 36 

Appendix B: Symbols 3 9 

Appendix C: Aerodynamic Data Employed in Investigation — 47 

Appendix D: Basis of Aerodynamic Model 77 



1 


INTRODUCTION 

The Langley Research Center has a broad research program designed to 
advance the state-of-the-art in the area of stall/ spin technology. One ma- 
jor requirement existing in this area is the development and validation of 
reliable theoretical methods for prediction and analysis of stall/ spin char- 
acteristics. Although theoretical studies have been made in the past, no 
concentrated nor continuing effort has been made to rigorously validate the 
analytical techniques employed in these studies. In view of the urgent need 
for valid theoretical techniques, this study was performed as a step toward 
advances in this area. 

The present investigation is the second phase of a two-part study di- 
rected toward determining the correlation between theoretically predicted 
and experimentally measured spinning motions of a-j^-scale, unpowered, 
radio-controlled, dynamically b all asted, drop model of a current fighter 
configuration. In the first phase, developed spin motions obtained with 
theoretical models based on conventional aerodynamics were compared with 
motions based on rotary (rotation balance) aerodynamics and factors, such 
as possible errors in aerodynamics and mass characteristics, control time 
histories, etc., which might influence the correlation between predicted and 
measured airplane spin motions were examined. This previous study reported 
in reference 1, which dealt primarily with the developed spin, presented the 
following conclusions, and consequently dictated the technical approach em- 
ployed in the present study: 

• The spinning motion can be acutely sensitive to a small change in 
the application of the spin entry controls. Since the control phas- 
ing required to trigger a radical change in the spin behavior of the 
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airplane cannot be predicted at this time, flight recorded con- 
trol time histories must be faithfully duplicated during an ex- 
perimental-analytical correlation study. 

• The spinning motion, both the incipient and developed phases, 
are sensitive to the location of the center of gravity but are 
not significantly affected by reasonable errors incurred in 
measuring airplane inertias. 

• The damping in pitch derivative, C , has a significant influence 

on all phases of the spinning motion. The magnitude of the static 
yawing moment in the 50 to 90 degree angle-of-attack range appre- 
ciably affects the incipient spin but only slightly affects the 
developed spin. 

• No spin can be computed using the conventional analytical technique 

when the dynamic (forced-oscillation) derivative C , C or C is 

m n 1 
q r * r 

unstable in the angle of attack region for spin equilibrium. When 
the contribution of the dynamic derivatives are limited to the 
oscillatory component of the total angular rates, a spin can be com- 
puted which, however, does not match the experimentally determined 
spin motion. 

• The spin computed with rotary data duplicated the developed spin 
obtained in the spin tunnel. Whereas, the spin computed in the con- 
ventional manner with static aerodynamic data was considerably less 
severe. It is indicated, therefore, that the aerodynamic moments 
generated in the spin due to steady rotational flow, as measured by 
a rotation-balance, are indeed significant. 
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• It appears that static aerodynamic derivatives can be extracted 
from rotation-balance data. It may be possible, therefore, to 
simply add the aerodynamics associated with steady rotation flow 
to the conventional aerodynamic model. 

The present investigation completed the two— phase study by determining 
the degree of correlation that can be obtained with measured flight motions 
by using a theoretical model based on the use of rotary, static and dynamic 
aerodynamics according to the technical approach identified in the first- 
phase study. To minimize Reynold's number effects, the experimental motion 
and measured aerodynamics were obtained with the same model. 
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TECHNICAL APPROACH 
Radio Control Model Time Histories 

The objective of the present study -was to correlate analytically pre- 
dicted model motions with the flight motions measured on a remotely-controlled, 
unpowered model launched from a helicopter and controlled from a ground-based 
pilot station. In this free-f light model technique, the model is launched 
and then controls are applied to force the model into the stall angle of 
attack range and possibly into a spinning motion. The test technique is used 
at the Langley Research Center for stall/spin testing. Parameters recorded 
during a typical flight include airspeed, angle of attack, angle of sideslip, 
the three body axis angular rates, and the control surface positions. Attitude 
angles are not currently measured. Using this model test procedure, the com- 
plete flight history of a vehicle entering a spin can be obtained, including 
the post-stall, the incipient spin, the developed spin, and the spin recovery 
motions . 

Ideally, for the purposes of a correlation study, one would like to have 
numerous flights wherein the pilot attempts to repeat the same spin-entry ma- 
neuver with the model to establish the repeatability of the model motions and 
to determine if the model exhibits any aerodynamic asymmetries which may make it 
more prone to spin in a particular direction. These repeated flights were not 
possible in the present study and therefore only one model flight was obtained 
which clearly exhibited the spin-entry, spinning, and spin-recovery motions of 
interest in this study. The time-history obtained during this flight is pre- 
sented in figure 1, and this is the flight which this study attempted to analyti- 
cally predict using the same control inputs that were applied to the flight model. 
The values of a, g, and V presented in figure 1 are referenced to the airplane 
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center of gravity and were computed by transferring the values measured at the 
model nose boom using the formulas in appendix A, where the values of AX, AY, 
and aZ were 5.17, 0.39, and 0. feet, respectively. 

Aerodynamic Model 

A large quantity of measured aerodynamic data, including static, dynamic, 
nnH rotary types, has been obtained on the subject model configuration at low 
speeds. A summary of the data available is presented in Table I showing the 
range of test variables covered. As can be seen, the data available were 
restricted in several variables. For instance, the rotation balance data was 
not measured below 0=55° and complete data were not available on each control 
surface for its full deflection range at high angles of attack. comparing 
the measured model flight excursions shown in figure 1 with the aerodynamic 
data available (shown in Table I), it can be seen that the aerodynamics were 
not measured to cover all of the angles of attack, sideslip angles and control 
deflections experienced in flight. Acknowledging these data limitations, an 
effort was made to construct the best aerodynamic model possible from the avail- 
able data. Probably the most severe limitation in aerodynamic data was with the 
rotation balance (or rotary) data which was restricted to only four control con- 
figurations and to angles of attack above 55°. 

The aerodynamic data referred to in Table I was extensively analyzed and 
combined, as discussed in appendix D, to develop the final aerodynamic model 
used for the analytical predictions. The aerodynamics used for the final model 
are presented in appendix C. These data represent the clean, scale model 
configuration used in the radio-controlled model flight test . All the data 
presented in appendices C and D were reduced to standard coefficient form on 
the basis of the following geometric characteristics: 


CABLE I: SUMMARY OF AVAILABLE AERODYNAMIC DATA FOR CONSTRUCTING AERODYNAMIC MODEL 
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S = 5.65 ft . 2 c = 11.76 in. b = 6. Ul ft. 

The data of Appendix C were transferred during the computation of time his- 
tories to a 20 percent c moment reference center to correspond with the flight 
model center-of-gravity location. 

The static and rotary aerodynam i c data measured on this particular model 
configuration indicate that aerodynamic derivatives such as 0^ * C , C , etc. 

e e i s 

can be obtained from the rotary data in a manner similar to that used with static 
data. Also, for this particular configuration, these derivatives and the longi- 
tudinal aerodynamics did not appear to vary greatly with rotation rate. Therefore, 
the conventional method of representing the longitudinal aerodynamics and the 
stability and control derivatives was used and the measured changes in the lat- 
eral-directional aerodynamics due to rotation rate were included in an incremental 
form. Although this approach seemed reasonably valid for the particular config- 
uration under study, it may not be valid for any airplane since rotary aerodyna- 
mics are known to be highly configuration dependent. However, if such a method 
can be used, as done here, it is very convenient since one can employ all of the 
static aerodynamic data available, and one can develop an aerodynamic model which 
can transition from the static aerodynamics valid at low rotation rates to the 
rotary aerodynamics needed at the high rotation rates and angles of attack 
associated with the spin entry and developed spin. Moreover, if the static con- 
trol data can be used, much more flexibility is available in simulating complex 
control time histories, and one is not restricted to the particular control con- 
figurations tested on the rotation balance. In view of the above discussion, 
the following aerodynamic model was selected for this study. 
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(o'.B ) 


(or»P ) 


(a»P jig) 


(a) 

(a) 

(a) 


Two aerodynamic models were employed herein which differed only in respect 
to the dynamic derivatives selected from Appendix C. One model employed the dy- 
namic derivatives that were experimentally obtained in the Langley full-scale 
tunnel facility at a nondimen sional forced-oscillation frequency of ^ = .0^7 
and the other used estimated dynamic derivatives. Unless noted otherwise, the 
experimental derivatives were always employed in this study. 

The data presented in Appendix C were inputed into the computer in tabular 
form. Values for control settings other than those presented were obtained by 

linear interpolation. By tabulating C , C and C as a function of <v and 8 for 

m n c H 

i =0 and -30°, the stabilizer control derivatives C , C, T and C were not 
s m. N. c. 

i x i 

s s s 

required. Interpolation of these tables would not incur errors except at low 

a where the tail is stalled for an i = - 30 ° setting-a flight condition not 

s 

to be simulated, during this investigation. If an angle of attack/or side- 
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slip angle were computed, momentarily -which was greater or less then those avail- 
able in the tables, the last values used would be maintained until the computed 
angle fell within the confines of the available data. 

The ability of this aerodynamic model to represent the radio control model 
requires that three assumptions be made. To wit: 

1) Rotation balance data does not change appreciably at lower angles of 
attack from that measured at 55° angle of attack. 

2) Control characteristics for this configuration are not influenced by 
rotational flow. 

3) The wing-mounted, upper surface spoiler used on the radio-control 
model had a negligible influence at high angles of attack. 

Inertia and Mass Data 

The following weight and inertia characteristics were used to represent 
the model : 

Weight = 145.12 lbs 

Center of gravity location = 0.20 c 

I = 1.805 slug-ft I 2 I =0. 

X xz 

I =7.326 slug-ft 2 

y 

I =9.388 slug-ft 2 

2 

These model values therefore represent a model which has relatively more mass in 
the fuselage than the wings, as indicated by the inertia parameters, i.e., 

I -I I -I I -I 

x ■ ■■£ =-.0298, =-.0111, = .0409. 

mu mo mb 
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Analytical Technique 

This study required, the computation of large angle motions using a 
large angle, six-degree-of- freedom, plotted output digital computer program 
with nonlinear, nrultifunctioned aerodynamic coefficients. This program solved 
the equations of motion and associated formulas presented in Appendix A as 
discussed in Reference 1 . 

Forty-second time histories presented herein were calculated employing 
the control time histories shown in Figure 1 for the radio-control model. 

Each time history is presented in a figure consisting of two 11 " x 17 " pages 
nnd each page contains the following variables plotted versus time. 

Page 1 Page 2 

a P 

8 1 

V R 

It was felt that the spinning motion was sufficiently illustrated by de- 
scribing the relative wind vector at the center of gravity on the first page and 
the an gul a/r velocities about the axes having their origin at the center of gravity 
on the second page. 

Unless noted otherwise, the following initial values at time equal zero 
seconds was employed for these studies. 


H' 

4250 ft 

P 

- 0.41 

rad/sec 

o 



V R 

110.2 ft/sec 

q 

0 

rad/ sec 

V 

-60.1 deg 

r 

0 

rad/sec 

oc 

26.0 deg 

0 

-50.8 

deg 

p 

0 deg 

9 

-50.0 

deg 

6 a 

0 deg 

Y 

0 

deg 

6 R 

0 deg 

i 

s 

-1.8 

deg 


I 
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All of these values except for the Euler angles -which are discussed later, 
were obtained from the flight records at a time equal 0 seconds as shown in 
Figure 1 which was the time at which the drop model was assumed to be out of 
the flow field of the helicopter. 
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RESULTS AND DISCUSSION 

The study of reference 1 had conclusively demonstrated the need for and 
the adequacy of the mathematical model used herein. The aerodynamic model was 
the best that could be developed under the constraints discussed in the Techni- 
cal Approach and Appendix D. Obviously, the level of confidence in the validity 
of the aerodynamic model was considerably less than in the mathematical model. 

The various parameters that were investigated, therefore, and discussed herein 
are identified in the List of Figures on page 19 . 

The ability to compute the flight recorded motion shown in Figure 1 when 
experimental and estimated dynamic (forced-oscillation) derivatives were used is 
shown in Figures 2 and 3, respectively. It can be seen that the experimental dynamic 
derivatives resulted in a considerably more oscillatory motion, which appears 
to diverge, than was evidenced by the radio-control model. When the estimated 
dynamic derivatives were employed, the computed motion was appreciably different 

in that the amplitudes of the oscillation were greatly reduced. In fact, the 
amplitudes in p, q, a and 0 were even less than was realized with the radio- 
control model. And, in this instance, the maximum magnitude of r was greater 
and the vehicle did not quite recover using the radio-control model recovery con- 
trol time histories . 

As shown in Appendix C, there is an appreciable difference between the 
experimental and estimated dynamic derivatives. However, since reference 1 had 
shown C to be of paramount importance and the magnitude for the experimental 

and estimated derivative was approximately -2 and -17, respectively above a = 75°, 

only the magnitude of the experimental C derivative was modified such that its 
" m 

value was -12 at and above = 75° ■> As shown in Figure 4, with this one change 
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incorporated in the experimentally determined aerodynamics model, the build- 
up in the oscillation shown in Figure 2 is eliminated and the oscillatory 
motion obtained between 12 and 24 seconds approximates the frequencies and 
amplitudes of the recorded radio-control model values. Consequently, the 
average values in B, p and q agree well, although the computed average a value 
is slightly less that the flight value (i.e., 87 and 91° for the computed and 
flight motions, respectively). Also, the flight recorded velocity and number 
of turns (turns not presented herein) were duplicated during the computed mo- 
tion. It should be noted that when the computed a trace oscillated to a peak 
value of 97°, only the aerodynamic values realized at a = 90 were employed 
throughout the time periods a exceeded 90°. This being the case since no 
aerodynamic data were available for angles of attack above 90°. Consequently, 
not matching the average oscillatory flight value of 91 is n °f surprising. 

Although the flight recorded spin and recovery motions were computed rea- 
sonably well, the incipient spin motion (i.e., from unstalled flight to a flat 
spin condition) that takes place in approximately the first 12 seconds after 
control movements are initiated does not correlate with the radio control time 
history. 

As mentioned previously, the flight attitude angles were unknown, lixfor- 
tunately, it was found that the initial motion was highly dependent on the 
assumed initial values for the attitude angles 9 and 0. For instance, large 
oscillations in a resulting in negative angles of attack for which no aerodyna- 
mic data was available and a full 360 ° roll could be computed during this time 
depending on the assumed initial values for 9 and 0. Since many combinations 
of these attitude angles are possible, the first combination which resulted in 
no violent a upset was used throughout the study even though they probably did 
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not represent the flight values the drop model had assumed at the time control 
movements were initiated. An error in the initial velocity can also produce 
this type of problem. 

Although it was known that the initial values of e and <t> could critically 
influence the resulting incipient spin motion, concern still existed relative 
to the differences between the flight recorded and computed motion. Consequent- 
ly, other factors were examined. For instance, knowing the importance of changes 
in control deflection, the stabilator deflection was increased from -12 to 20° 
between 6.80 and 12.18 seconds. The motion computed with this control change 
is presented in Figure 5. A significant delay in the time for r to increase 
was realized such that the magnitude of r attained at the time of applying re- 
covery controls was reduced from 13.4 to 11.6 rad/sec. Although this value is 
closer to the flight recorded value of 12 rad/sec, the overall r trace correla- 
tion is further degraded. 

It was also thought that the computed r trace during the incipient spin 
phase was the consequence of applying arbitrarily the rotational flow charac- 
teristics obtained at <* = 55° to lower angles of attack. The effect of limiting 
the influence of rotational flow to angles of attack above 55° was, therefore, 
investigated. It was found that in this instance the incipient spin phase 
correlation deteriorated as did the recovery phase (compare Fig. 1 with Fig. 

6). It was indicated, therefore, that rotational flow effects are important 
below a = 55° and that in the future rotational flow data should be obtained 
down to at least 30 ° angle of attack. 

The control time history employed on the radio control model did not trig- 
ger a spin in the shortest time possible nor was the flat spin completely sta- 
bilized (relative to yaw rate) before recovery controls were applied. This 


i 
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situation was realized while examining the effect of changing the stabilator 
deflection between 6.80 and 12.18 seconds. It was observed that when the 
stabilator deflection was increased from -12 to -25 during this time period, 
the incipient spin motion did not culminate in a flat spin to the right as 
the model did. However, when the recovery controls were introduced beginning 
at 17.9 secs, (which previously resulted in a recovery from a right spin), the 
vehicle went very rapidly from a post-stall gyration motion to a spin to the 
left. The resulting motion which ensues when the recovery time his- 
tory is introduced at t = 0 seconds is shown in Figure 7* Also, the initial 

values in this instance were 


V R 

l6l.8 ft/sec 

P 

- 1.57 

rad/ sec 

7 

-48.5 leg 

q 

0 

rad/sec 

CC 

0 deg 

r 

0 

rad/ sec 

0 

8.1 deg 

0 

- 82.5 

deg 



e 

- 56.5 

deg 


which were the values existing at 17-9 seconds in the above described motion. 

It will be noted that some of these values are different from those previously 
prnp l oypfj . As can be seen in Figure 7 , the duration of the incipient spin phase 
has been considerably reduced and, as shown, the shape of the r trace is more 
closely simulated in this instance. This may be the result of assuming dif- 
ferent initial values for 0, 0 and V and/or the fact that very little time 
is spent in this instance below « = 55 0 for which rotational data had to be 


assumed. 
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CONCLUDING REMARKS 

This correlation study was conducted to determine what degree of corre- 
lation could be obtained between analytically predicted stall/spin motions and 
motions measured during flight tests of an unpowered, remotely-controlled scale 
model of a current fighter configuration. The emphasis in the study was to 
determine if the use of rotation balance aerodynamic data would provide improved 
correlation. Due to time limitations, only one model flight time history was 
available for the correlation study. This fact prevented this study from de- 
termining the degree of repeatability of flight motions that could be obtained 
in the model flights. The lack of this information prevents an evaluation of 
possible errors in correlation due to such factors as model aerodynamic asym- 
metries and small variations in control deflection timing. 

The results of the correlation study showed that the aerodynamic model 
formulated to use the rotation balance aerodynamics was necessary to obtain the 
best possible prediction of the model flight motions. Specifically, the corre- 
lation was best for the developed spin motion and relatively poor for the spin 
entry and final phase of the recovery motion. In this regard, it would appear 
desirable to measure rotational flow effects down to lower angles of attack 
than is presently done, at least down to near a = 30°. Furthermore, more work 
is needed to refine the aerodynamic model needed to represent the motions 
occurring during the spin-entry motion. Several other important factors identi- 
fied in the correlation study included: 

(l) A dominant effect of the value of the damping parameter 
C in the developed spin region. 

m q 


i 
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(2) A major influence of the initial values of the Euler angles used 
at the initiation of the motion calculation - these results indi- 
cate that it is necessary to accurately know the initial flight 
attitude of the vehicle when trying to analytically predict the mo- 
tions. Therefore, the attitude angles of the flight vehicle must he 
accurately known. 

In summary, the major result of this correlation study is that rotation balance 
aerodynamic data are needed to formulate a useful analytical model for spin 
calculations . Shortcomings of previous analytical correlation efforts, which 
attempted to use the conventional static and dynamic stability derivatives, 
were probably ultimately related to the omission of rotary aerodynamics. A 
final caution to the reader regarding the application of these results to other 
vehicle configurations is that it is well-known that the rotation balance aero- 
dynamics measured to date on a range of airplane configurations show the rotary 
aerodynamics to be quite configuration dependent. For instance, configurations 
other than the one studied herein may very well exhibit variations in both static 
stability and control power with rate of rotation. Therefore, due caution should 
be taken in extending the assumptions of this study to different configurations 
without careful study of the particular data. 
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APPENDIX A 

EQUATIONS OF MOTION AND ASSOCIATED FORMULAS 


The dy n a m ical equations required to specify the translational and rota- 
tional motions of a rigid body moving through space are described in this 
appendix. The familiar six degree of freedom differential equations repre- 
senting linear and angular accelerations of a moving body axis system having 
its origin at the aircraft center of mass are presented below. 


E F 

xl = - g sin 0 + vr - wq + aero 

m 


v = g cos G e sin d + wp - ur + aero 

m 


EF 

w = g cos 0 cos <f> + uq - vp + aero 

e e 

m 


P = ~*y ~ qr + X xz (r + pq) + ^aero 


EL 


I-I I ,2 2s , jM 

- xz (p - r ) + aero 


q = z x pr 

I I 

y y 


I - I , I fj, s EN 

- + xz (p - qr) + aero 


r = x y pq 

I I 

z z 


In addition, the following formulas were used: 


a - 

tan 

- 1 ® 



u 


B = 

sin' 

•1 / V 

\ 

l‘v~ 

1 



\ R 

/ 

V = 

jf~2 

2 ^ 

2 

R 

U u 

+ V + 

w 

XL - 


^ 2 4 

2 

Mp 

+ q + 

r 


i 
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Turns in spin = 


JV* 


* =. ° > 
T e sin 0 


( t .-If sin 
e = sin l — — 
\ cos 


fi = q cos <j> - r sin 0 


0 = p + r tan 0 cos 0 + q tan 9 sin 0 

6 6 e e 


P = p + P^ 
r o 


<1 - 9* + * c 


r = r + r 
r o 


These total angular velocities (p, q, r) consist of steady rotation (p^, q^ , 
r ) components upon which oscillatory (p Q , q^, and r^) components are superimposed. 
These components are defined as follows : 


p = -7 sin 0 

r e e 

q = t cos 0 sin 0 

’r e e e 

r = ^ cos 0 cos 0 

r e e e 


p = 0 

o e 

q = 0 cos 0 

o e e 

r = - 0 sin 0 
o e e 


For the aerodynamic model, the following total derivatives were used: 


C 1 = C + C V 

N w \-w- 

c 1 = c 

c c 

c 1 = c + c 8 + c a + c + C V? + C V 3 

y y y 6 a y 6 R y rot y r 2V y p 2 V 

a R 

C 1 = C+C 6+C 6+C +C r o b + c P o b 

* a ^6 a R *rot A r 2V *p 2V 

a R 
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Cl 

n 


C + C 
n 


n. 


6 + C 


6„ + C 


+ C 


n. 


R 


rot 


n. 


+ C 


n 


P o* 

2V 


C* 

m 


C + C 
m 


m 


q c 

2o 

2V 


Note: 


The values of C 

y rot 

using the values of 


, C and C are obtained from rotary balance data 
*rot n rot 

q^ and r_^ to compute ^ as 



The following formulas were used to compute a and g at the nose boon vane 
location. 


^oom 

v 

Boom 

W Boom 


u + Az.q - Ay .r 
v + Ax*r - Az.p 
w - 12 c.q + Ay*P 


tJ 2 2 ,2 

V7’ = u + v + 

Boom ±$oom Boom 

“Boom = tan _1 -£222. 

^Boom 


Bt, = sin _1 v ] 
D Boom — 


r 

Boom 


2 

w 

Boom 


where Az , Ax and Ay are the distance between center of* gravity and boom vane 
location. 
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APPENDIX B 


SYMBOLS 


4o 


Measurement s and calculations were made in the U.S. Customary Units. Fac- 
tors for converting these units to the International Systems of Units may be found 
in reference 2 . 

The Body Axes System is used with the origin of the axes system located at 
the aircraft center of gravity. The X axis is parallel to the aircraft Fuselage 
Reference Line (FRL) and is positive forward, the Y axis is positive towards the 
right wing tip and the Z axis is positive downward. The following illustration 
shows this axis system and the positive direction of the angles, forces and 








In the following definitions a dot ( . ) over a quantity denotes one 
differentiation with respect to time and a zero subscript (o) denotes the 
initial value of the quantity. 


Xj, Y , Zj , inertial axes 
Xg, Y b , Zg "body axes 



BODY AXIS SYSTEM LOCATED IN INERTIAL SPACE 
POSITIVE DIRECTION OF EULER ANGLES AND RATES 
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I 

z 


I 

xz 


m 

W 

S 


b 

c 

W/S 

dy - i z )M 2 

(I 2 - I x )/»b 2 

(I - I )/mb 2 
x y y/ 

t 


P 

q 

H 


0 


V 




e 


0 


p 

moments of inertia about the X,Y and Z slug-ft 

body axes, respectively 

O 

product of inertia, positive when the slug-ft 

principle X axis is inclined below the 
X body axis at the aircraft nose 


aircraft mass (=W/g) 
aircraft weight (=mg) 
wing area 
wing span 

wing mean aerodynamic chord 
wing loading 

inertia rolling moment parameter 
inertia pitching moment parameter 
inertia yawing moment parameter 
time 

atmospheric density 

V 2 

dynamic pressure = 
vertical height 
angle of attack 
angle of sideslip 
Euler angles 

angle between Y body axis and horizontal mea- 
sured in vertical plane, positive for erect 
spins when right wing is downward and for 
inverted spins when left wing is downward 


slug's 

lbs 

ft 2 

ft 

ft 

lb/ft 2 


sec 

slug/ ft 

lb /ft 2 

ft 

deg 

deg 

deg 

deg 


i 
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Y 


i 

s 


6 


R 


V. 


R 


VL, V, W 


U, V, w 


p, q> r 


0) 


XL 


flight path angle, measured, in a vertical deg 

plane between the horizontal plane and the 
resultant velocity vector, positive when the 
resultant velocity vector points above the 
horizontal plane 

differential horizontal tail deflection (half deg 
angle), positive to produce left rolling mo- 
ment 

control deflection-stabilizer, positive direc- deg 
tion trailing edge down 


control deflection-rudder, positive direction deg 
trailing edge left 


total linear velocity vector 
free- stream velocity 

components of the total linear velocity 
vector (V ) along the X, Y and Z body 
axes , respectively 
components of the relative linear 
acceleration vector along the X, Y 
and Z body axes, respectively 


ft/sec 

ft/sec 

ft/sec 


^ / 2 
ft/sec 


components of the total angular velocity 
vector along the X, Y and Z body axes, 
respectively 

frequency of forced-oscillation tests 
total angular velocity vector 


radians/ sec 


radians/sec 
radians/ sec 
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g 


F ,F , 

X V 

aero J aero 


aero 

L , M , 
aero’ aero’ 


N 


aero 


-F 

x 

1 _ aero 

V ' 

c q S 


1 ^aero 

u — 1 
y of? 


-F 


c 1 = 

N 


aero 


q S 


p ± _ aero 

U 1 " — 

qSb 


M 


m 


aero 


qSe 


1 N 

q x _ aero 
n 

qSb 


C c ^ a * i s^ 


C y (P, a) 


2 

acceleration due to gravity ft/sec 

components of the aerodynamic forces along lbs 

the X, Y and Z body axes, respectively 

components of the aerodynamic moments ft 'lbs 

about the X, Y and Z body axes, respectively 

total engine thrust force lbs 

total chordwise force coefficient, 
positive along the negative X body axis 

total sideforce coefficient, positive along the positive 
Y body axis 

total normal force coefficient, positive along the 
negative Z body axis 

total rolling moment coefficient, positive direction 
of moment drives the right wing tip down 

total pitching moment coefficient about the aerodyna- 
mic center, positive direction of moment drives the nose 
up 

total yawing moment coefficient, positive direction of 
moment drives the nose right 

variation in chordwise force coefficient due to angle of 
attack, sideslip angle and stabilizer deflection 

variation in side force coefficient due to sideslip 
angle and angle of attack 


i 
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C (/Lb/2V,a,i ) 
y rot S 


N 


(a» 8, i g ) 


variation in side force coefficient due to rate of 
rotation, angle of attack and stabilizer deflection 

variation in normal force coefficient due to angle of 
attack, sideslip angle and stabilizer deflection 


C g (fl, a, i ) 

Z s 


C Ub/2V, 
£ + 



C (a, P, O 
m s 


c (p, a, i ) 

n ^ s 


variation in rolling moment coefficient due to side- 
slip angle, angle of attack and stabilizer deflection 

variation in rolling moment coefficient due to rate of 
rotation, angle of attack, and stabilizer deflection 

variation in pitching moment coefficient due to angle 
of attack, sideslip angle and stabilizer deflection 

variation in yawing moment coefficient due to side- 
slip angle, angle of attack and stabilizer deflection 


C 

4 - 

rot 


(j\b '2V, 




(“» i s ) 


6 


i 

s 


) 


c 

n 


6 




a 


i ) 

s 


(oh 

6 


l 

s 


) 


variation in yawing moment coefficient due to rate of 
rotation, angle of attack, and stabilizer deflection 

variation in due to angle of attack and stabilizer 
* 6 a 

deflection, per deg 

variation in ^ y due to angle of attack and stabilizer 

*iT 

a 

deflection, per deg 

variation in ^^n due to angle of attack and stabilizer 



deflection, per deg 

variation in ^ C y due to angle of attack and stabilizer de- 
flection , per deg 
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C 





ck? P ) 


n. 


(i 


Of> 


0 ) 


C 


Z 


N 




9 


£b 

2V r 



c 



= 9 C 

IE 

2V r 
= 9 C 

IE 

2V r 




r*C 

variation in fl due to stabilizer deflection, angle of 

96 r 

attack and sideslip angle, per deg 
9C 

variation in n due to stabilizer deflection, angle of 

*«R 

attack and sideslip angle, per deg 
per deg 
per deg 

per radian 
per radian 
per radian 
per radian 

per radian 

per radian 
per radian 

per radian 
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Figure C-9. Effect of angle of attack and sideslip angle on yawing-moment coefficient 
when i„ = 0° and -30°. 
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APPENDIX D 

BASIS OF AERODYNAMIC MODEL 
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A considerable amount of effort was expended in defining the aerodynamic 
model presented in Appendix C. This effort basically involved the following 
three tasks: 

1) Analyze the Langley full-scale tunnel data to identify the significant 
parameters affecting the aerodynamic characteristics and extract the most valid 
model from these data. 

2) Validate the use of Langley 12-foot data for modeling control characteristics not 
directly available from Langley full-scale tunnel tests. 

3 ) Identify model required to represent rotational rate effects. 

Langley Full-Scale Tunnel Data 

Figures D-l through D-7 identify the influence that angle of attack, sideslip 
angle and control configuration may have on the longitudinal, lateral and direc- 
tional stability characteristics. It can be seen from these figures that the 
stability characteristics are significantly effected by & , fi and stabilizer 
deflection and not by the rudder or differential tail setting. The variation 

of C, C and C with oc, and /3 presented in Appendix C was obtained, therefore, 
i ’ n y 

by fairing through the data obtained for the different control configurations 

associated with a i =0 or -30 setting. 

s 

The curves of Figures D-L through D-7 were obtained from Figures D -8 through 
D-25. These later figures present rolling- and yawing-moment and side-force 
versus cc for various yfl angles and control configurations. In all instances, 
the sign of the force or moment was reversed for positive ft values. (The 
closeness of the corresponding positive and negative JS data for lateral controls 
neutral, as well as the data points for Q = 0 ° themselves, demonstrate model 
symmetry.) Curves were then faired between the plotted positive and negative 
data point. These curves were then assembled and compared in Figures D-4 

through D-7. Fairing the data in this manner allows one to eliminate the control 
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moment and force increments and to examine the effects of control deflection 
on stability. Therefore, the absolute increments due to control deflections do not 
appear in Figures D-4 to D-7. 

Rudder and differential tail control characteristics at sideslip angles 

were also determined from these figures by taking at a given c£ one half the 

incremental value between corresponding positive and negative sideslip data 

points from the appropriate control configuration figure. (The value for fi =0° 
could of course be read directly off the figure.) These control characteristics 

were then compared with the values obtained by substracting the neutral controls 

configuration data from the subject control configuration data. 

Langley 12-Foot Low Speed Tunnel Data 

Given below are the control configurations for which lateral-directional 
data were available from the Langley full-scale (30 , x6o’) and 12' facilities. 
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*Also rotary balance data available from = 55 to 90° • 


As can be seen, the isolated effect of rudder and differential tail deflection at 
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high values of stabilizer setting were not available from the 30 , x60* facility. 
Consequently 5 it was desired to employ the available 12 * control configuration 
data in this area. 

The curves presented in Figures 26 and 27 were prepared using the same tech- 
nique described above for the Langley full scale (30 l x60 1 ) data analysis. It 
can be seen from these figures that, as was the case for the 30'x60 T data, 
the rudder and differential tail did not effect the lateral-directional sta- 
bility characteristics. However, the rolling and yawing moment characteristics 
obtained in the 30'x.60' and 12 1 facilities for the identical i g = 0° control 
configurations (presented in Figures D-28 and 29, respectively) do not correlate 
well over some angle of attack regions. That is, less lateral stability in the 
& = 25 to 40° region and less directional stability above Of = 30 ? especially 
at high JH values, is realized in the LRC 12' facility. Obviously, some model- 
tunnel wall interference effects were encountered with the Jq ~ scale model. 

It was decided, therefore, not to use the 12 1 data directly but to obtain in- 
cremental effects due to control deflection and to determine the influence of 
stabilizer setting on other control characteristics. It was found that rudder 
control characteristics obtained from both facilities did indeed compare favor- 
ably (see Figure D-30) and it was shown that the i fi setting did not influence 
rudder power (see Figure D-3l) . The differential tail control power from both 

facilities was also in good agreement but i settings (as shown in Appendix C) 

s 

had. an appreciable influence on the control power. 

Rotary Balance Data 

Rotary balance data presented in reference 1 for clockwise and counter 

clockwise rotation -O- b values of .04-93, .1030, .1523, *3015, .2553 and .30^5 

2V o o 

at d and (5 values of 55, 60, 65, 70, 75, 80, 85 and 90 and 0 and 5 , re- 
spectively, were analyzed for the following available control configurations : 


8 o 


is’ »a’ 6 R 

= 0° 




i s - -30°, 

- 

0°, 

6 = 

R 

0° 

i s = -30°, 

4 a = 

5°, 

*R " 

-30° 

i s = -3°°> 

«. * 

-5°, 

O 

Sd 

ll 

30° 

By connecting the clockwise 

and 

counter 

clockwise data points with a straight 


line segment, forces and moments at^ ± 0 were determined. The static longitudinal 
characteristics obtained from the rotary balance data using this procedure are pre- 
sented in Figure D-32. The i g = 0 and -30 control configuration data of Figure D-32 
are compared in Figure D-33 with static test data. It can be seen that although the 
C m’ C N an<i C c varia ' bions with of are approximately the same (also illustrated in 
Figure D-34 which compares the control characteristics), they are displaced by 
approximately a constant value. The magnitude of C m and C N being less and C c greater 
for the coefficients extracted from rotary balance data when compared with the 
static test data. Since the longitudinal characteristics were invariant for this 
configuration with rate of rotation, these incremental effects would have been ex- 
perienced between no rotation and an || value of .0^93 (first test point). This 
did not seem probable, therefore, these incremental effects were not included in 
the aerodynamic model. It was found that the lateral and directional stability 
characteristics were in fair agreement (see Figure D-35) with the values established 
during the static tests. 

The longitudinal characteristics for this model were found to be essentially 
invariant with rate of rotation and, therefore, it was not necessary to include any 
effect of rotation rate in the longitudinal aerodynamic model. This was not the 
case relative to C^, and C^. To establish the influence of rotation rate, the y5 
= 0 counter clockwise data vs ~ were fitted over the clockwise data points and 
the same procedure was applied separately to the/3 . 5 0 data. Then the - pj^data for 
A= 5° were adjusted by the incremental value obtained between/? = 0 and 5° at ~ = 0 


I 
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and superinposed on the y3 = 0° data. Therefore, four sets of data points vs 

were available for each angle of attack and control configuration through which 

slopes could be drawn. It was found that all the i g = -30° control configurations 

has the same rotation rate characteristics which, however, were different from the 

i = 0^ characteristics. Consequently, rotation flow models are presented in 
s 

Appendix C for both i = 0 and -30 settings . 
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Effect of sideslip angle on pitching-moment coefficient for 
Figure D-3. ■ varlous control configurations. 





is =-30°, 6 a = 5°, 6 r 
F igure D-3. Concluded 





Figure D-4. Concluded. 






a) 0=10°, 20° *■ 

figure D-5. Effect of i g = -30° control configurations on yawing moment 
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-30° configurations 
Figure D-7. Concluded. 
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Figure D-8. Effect of angle of attack and sideslip angle on yawing- 
"nnpnt, coefficient for configuration i st 6 n , 6 r = 0°. 
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Figure D~l4. Effect of angle of attack and sideslip angle on rolling- moment 

coefficient for configuration i s , 6 a , 6n = 0°. 
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Figure E-EO. Effect of angle of attack and sideslip angle on side-force 

coefficient for configuration i_, 6 a , 6u = 0°. 
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Figure D-21. Effect of angle of attack and sideslip angle on side-force 

coefficient for configuration i s = - 30 ; 6 a , 6 r = 0 . 
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Figure D-22. Effect of angle of attack and sideslip angle on side-force 

coefficient for configuration i s , 6 r = 0°; 6 a = 12°. 
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Figure D-H4. Effect of angle of attack and sideslip angle on side-force 
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